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One of the advances in computational design has been the development of adjoint meth- 
ods allowing efficient calculation of sensitivities in gradient-based shape optimization. This 
paper discusses two new applications of adjoint methodology that have been developed to 
aid in sonic boom mitigation exercises. In the first, equivalent area targets are generated 
using adjoint sensitivities of selected boom metrics. These targets may then be used to 
drive the vehicle shape during optimization. The second application is the computation 
of adjoint sensitivities of boom metrics on the ground with respect to parameters such as 
flight conditions, propagation sampling rate, and selected inputs to the propagation algo- 
rithms. These sensitivities enable the designer to make more informed selections of flight 
conditions at which the chosen cost functionals are less sensitive. 


Nomenclature 


p 1 + V 

/3 n Lagrange multiplier for absorption 

T Non-dimensional thermo-viscous parameter 

7 Ratio of specific heats, 1.4 

7 o,ru 7 i,n Lagrange multipliers for relaxation 
A n Lagrange multiplier for nonlinearity 

p Atmospheric density 

a Non-dimensional distance 

r Non-dimensional time 

r Retarded time 

0 u p-2 Non-dimensional relaxation time parameters 

A n ,B n Matrices during the first relaxation process 
A^ 2 , Blf Matrices during the second relaxation process 
A 3 , L >3 Matrices during the absorption process 
A e ,cp Equivalent area spline control points 

c Speed of sound, m/s 

a,i _ 2 Non-dimensional dispersion 

D Vector of design variables 

F F-function 

G Ray tube area 

h Cruise altitude 

In Objective function for adjoint calculation 

k n Scaling factor due to ray-tube spreading and stratification 

L Lagrangian 

M Cruise Mach number 

N Number of steps during propagation 

p, P Pressure waveform during propagation 

g, r, t Intermediate pressure waveforms 
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S Number of points used in propagation 

Xcp Axial spline control points 

Subscripts 

0 Starting conditions 

n Propagation iteration counter 

I. Introduction and Motivation 

There has been a lot of interest for a viable commercial supersonic aircraft after the successful demonstra- 
tion of the DARPA Shaped Sonic Boom Demonstrator (SSBD) 1 program. To mitigate the adverse effects 
of sonic boom, several efforts 2,3,4 ’ 5 , 6, 7 p ave p een carried out to design aircraft that overcome stringent 
acceptability constraints that currently forbid overland supersonic flight. Significant advances in computa- 
tional design in terms of tools and methodologies have been developed that complement the existing body 
of knowledge from the past several decades. One of the first steps in this design process is developing a 
framework that allows designers to explore the design space efficiently and develop concepts using appropri- 
ate analysis fidelity. In addition to non- gradient optimization methods, computational approaches based on 
sensitivities offer an elegant way to reduce sonic boom levels. Adjoint-based methods in particular provide 
an efficient way to compute sensitivities of chosen metrics to a large number of shape design parameters. 
Previous studies have demonstrated the capability of adjoint-based methods to optimize near-field pressure 
waveforms, 8,2,9 ground-based boom metrics, 10 as well as equivalent area metrics. 5,11 

The methodology described in this study builds upon the adjoint methodology from previous papers 10,11 
and extends the adjoint formulation in sBOOM 12 it in two ways. The first extension allows generation 
of equivalent area targets from sensitivity information which enables designers to quickly develop smooth 
targets for under-track and off-track azimuths in the neighborhood of the baseline design. The second 
extension calculates adjoint sensitivities of boom acceptability metrics to flight conditions (e.g., Mach, cruise 
altitude), propagation parameters (absorption and relaxation constants, step-size), and sampling rate used 
during propagation. This work facilitates uncertainty analysis and robust design of low-boom aircraft. 

II. Extension to Existing Theory 

The underlying mathematics behind boom adjoints is similar to that presented in Refs. 10 and 11. In 
this paper, we focus only on specific details that pertain to target equivalent area generation and sensitivities 
with respect to atmospheric and propagation parameters. 

II. A. Target equivalent area generation 

The theory behind computation of sensitivities needed in this case is a simple extension of the theory 
presented in Ref. 10. The primary difference is that the independent design variable vector is the input 
equivalent area instead of the off-body pressure waveform. To extend the sensitivities to equivalent area, 
we employ the chain rule of differentiation as shown in Eq. 1, where the relation dp/p = ^j==F is used. 
The first term on the right hand side is calculated from the previous formulation. The F-function (F) 13 
is calculated from equivalent area ( A e ) as given in Eq. 2. Numerically differentiating this relation yields 
which can be used in Eq. 1 to obtain the relevant sensitivity of the cost functional with respect to the 
equivalent areas. 


dL _ dL d(dpp) _ dL 7 M 2 * dF 
dAe d(dpp) dAe d(dpp) a/2 /3R dAe 

1 C y Ap" 

(2> 

The objective chosen to minimize is the A-weighted loudness of the ground signature. While it is desirable 
to use the perceived loudness level as the objective during optimization, there was no easy way to compute 
the sensitivity of this metric with respect to the ground signature that is needed in adjoint-based design 
optimization. A-weighted loudness computation, on the other hand, offers an analytical calculation procedure 
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that can be differentiated for generating required sensitivities and therefore amenable to adjoint-based design. 
In addition, the A-weighted loudness is well correlated 14 with the perceived loudness values. If a technique 
were to become available that would allow computation of perceived loudness level sensitivities, then the 
current methodology could be extended to minimize a cost functional that is a function of the perceived 
loudness of the sonic boom ground signature. 

At this point, a gradient-based optimizer can be used to change the baseline equivalent area values 
at different longitudinal sections along the aircraft to minimize the A-weighted loudness. However, this 
optimization can lead to equivalent area shapes that are not smooth. This is because the optimizer exploits 
shock cancellation in a favorable way due to the nature of the propagation algorithm. This is demonstrated 
in Figures 1, 2 and 3. Figure 1 depicts an optimized equivalent area that is visually smooth, but has ripples 
throughout when closely inspected. Looking at the equivalent area in F-function space reveals that the 
equivalent area ripples manifest as large oscillations around a smooth profile. During propagation, these 
oscillations cancel one another and the resulting ground signature improves relative to the baseline (See Fig. 
3) as predicted by the optimizer. 



Figure 1: An example of non-smooth target generated without spline smoothing. 

While non-smooth targets are numerically better, designers seek smooth targets not only for better 
performance but also for improved robustness in the final designs. To overcome non-smoothness in the 
target equivalent areas, an additional step is needed wherein smooth equivalent area shapes are generated 
by manupulating the control points of a spline. A similar approach using Bezier curves defined by arbitrary 
number of control points to represent the target equivalent area distributions is employed in Ordaz et al. 15 
One difference here is that the spline-based algorithm not only returns the spline interpolation but also 
produces the Jacobian matrices qxcp anc ^ dA^cp ' These matrices are then used in Eq. 3 to produce the 
sensitivity of the objective functions to the spline control points. In this paper, we generate a cubic spline 
and the corresponding Jacobian matrices. However, the formulation can be extended to other curve fitting 
approaches. 


dL 

dD 


dL dA e 
dA e dXcp 
dL dA e 
dA e dA e ,cp 


if D = X C p 
if D = A e) cp 


(3) 
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Ap (psf) F-function 



Figure 2: An example of non-smooth F-function generated without spline smoothing. 



Figure 3: An example demonstrating numerically better ground signature without spline smoothing. 
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II. B. Sensitivities with respect to flight conditions and propagation parameters 

The discrete adjoint equations derived in this section are based on a similar implementations developed in 
Ref. 10. The primary difference is the independent variable vector D, which has been modified to represent 
the 9-tuple vector [Act, S', T, C^i, h, M]. The formulation is not limited to this vector; other 

variables such as temperature and wind profiles may also be included. However, in this study we are limiting 
ourselves to the nine variables listed. 

The Lagrangian corresponding to a chosen objective may be written as in Eq. 4, where k n is the scaling 
factor given by • The details of the matrix terms are provided in the Appendix. Taking the derivative 

of the Lagrangian with respect to D results in Eq. 5. Contrary to the previous formulations, the relevant 
matrices during propagation now depend on the design variable vector D and hence the matrix derivative 
terms now appear in the derivative of theLagrangian. 


N 


N 


L(p, q, r, t, D) = I N (p N , D) + ^ 7 ^ n [ A n q n - k n B n p n 3 + ^ 77 ™ [A^r n - B n q n \ 


n= 1 


N 


n= 1 
N 


( 4 ) 
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dD dpN dD 

' / y ^0 ,n 

n=l 


8A- dp,,., dk , 

<j[) <)T> n dD dD^ P n dD Pn 


N 
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dD dD 


dD dD 


Qn 


N 


n = 1 


’ A n^n + d^_ tn _ B n dr n dB% 


dD dD 


dD dD 


( 5 ) 


N 


n= 1 


dp n df? dt n 


dD dtn dD 


The adjoint equations and their solution procedure are identical to previous formulations and are not 
repeated here. Once the adjoint terms are removed, the adjoint gradient can be calculated using the ex- 
pression in Eq. 6. It should be noted that the matrix Jacobian terms are three-dimensional matrices. It 
is clear by looking at this equation that the memory requirement increases significantly from the previous 
formulation. This is due to multiple matrix Jacobians and intermediate pressure waveforms that need to be 
stored at each step during the propagation process in addition to the Lagrange multipliers. Sensitivity of k n 
and consequently the sensitivity of the Blokhintzev invariant (See Eq. 7, in the absence of winds) 16,17 also 
needs to computed and stored. 


dL 

dD 


N 

J2d0,n 

n = 1 


'dA n ~ dk n „ n ~ u dB n _ 
dD Qn dD B Pn ~ l n dD Pn_1 


N 

n = 1 


dA$ 

dD 


r n ~ 


dD 


q n 


N 


Y.A 

n = 1 


93 

dD 


t n 


dB% 

~dD 


(6) 


P 



constant 


( 7 ) 


The matrix Jacobians can be obtained by differentiating the individual terms of the matrices in the 
Appendix. For sake of brevity, only a few of these are shown in the Appendix. 
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III. Verification of Adjoint Sensitivities 


To verify that the sensitivity values obtained from adjoint calculations are correct, a complex variable 
version of the process is developed. The use of complex step approach 18, 19 is standard practice 20 for gradient 
verification. The relevant algorithms are modified to work with complex variables and the derivatives of 
appropriate cost functionals with respect to the chosen design variables are calculated using an imaginary 
step size of 10 -50 . 


III. A. Target equivalent area generation 


Table 1 shows a comparison of the 
A-weighted loudness sensitivity to 
equivalent area values at selected 
indices, computed using the adjoint 
and complex versions of the code. 
It can be seen that the values agree 
to 9 decimal places, verifying that 
the sensitivity values are accurately 
computed. 


Table 1: Comparison of the adjoint and complex variable gradients 
for A e sensitivity. 


Grid Point 

Adjoint Gradient 

Complex Gradient 

0 

4.99555122854739 

4.99555122485994 

200 

0.22333517023795 

0.223335170026696 

500 

-1.83336776546880 

-1.83336776671720 


III.B. Sensitivities with re- 
spect to flight conditions and 
propagation parameters 


Table 2 compares the adjoint 
and complex variable gradients of 
ground loudness with respect to 
each of the flight condition and 
propagation parameters. The re- 
sults obtained using the adjoint im- 
plementation exhibit good agree- 
ment with the complex- variable ap- 
proach. Contrary to previous for- 
mulations, this sensitivity calcula- 
tion involves summation over the 
entire propagation distance, caus- 
ing some precision loss in the sen- 
sitivity comparisons. Nevertheless, 
the adjoint sensitivities are reli- 
able for any subsequent analysis. 


Table 2: Comparison of the adjoint and complex variable gradients 
for flight and propagation parameters. 


Variable 

Adjoint Gradient 

Complex Gradient 

s 

-0.121609572658E-002 

-0.121609679831E-002 

Act 

-0.101433622238E+004 

-0.101434351727E+004 

r 

-0.242603369041E+005 

-0.242607576473E+005 

c Vl 

0.850297505015E+002 

0.850293279821E+002 

C V2 

0.114275582996E+003 

0.114279749352E+003 

Out 

-0.270434382922E+004 

-0.270438763849E+004 

6^2 

-0.226512310051E+005 

-0.226534988755E+005 

M 

0. 103286198407E+002 

0.103281674303E+002 

h 

0.373155218390E-004 

0.373155285384E-004 


IV. Results 

IV. A. Generation of target equivalent area 

In this section, a gradient-based optimization process of generating smooth equivalent area targets is demon- 
strated. A baseline concept is evaluated at multiple azimuthal angles and its off-body pressure waveform 
and the corresponding area distributions obtained. Using the end areas, a cubic spline is fit to the equivalent 
areas at each of the azimuthal angles. In this demonstration, 15 control points are used for each spline 
fit. An SQP algorithm 21 was used as the gradient-based optimizer. Figure 4 depicts the design and target 
equivalent areas for each of the azimuthal angles. Smooth targets are generated for each azimuthal angle in 
the neighborhood of the underlying design area. 

A typical convergence history of the optimizer is shown in Fig. 5. Convergence is achieved in 50-60 
function evaluations, taking roughly 30 minutes of wall-clock time on a 64-bit Intel Xeon processor with 
16GB of memory. Targets at different azimuthal angles can be obtained simultaneously allowing designers 
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Figure 4: Generation of targets via gradient-based optimization using adjoint sensitivities. 


to quickly generate equivalent area targets. In addition, because the targets are in the neighborhood of the 
underlying areas at each of the azimuthal angles, there is a high degree of correlation between the targets. 
This can be helpful during shape optimization to make the problem better posed for obtaining an outer mold 
line that can match the under- and off-track targets simultaneously. Because of the fast turn around time, 
targets can be updated as the design evolves. The ground signatures corresponding to the generated targets 
are shown in Fig. 6. 

IV.B. Sensitivity to flight conditions and propagation parameters 

Figure 7 shows the sensitivity of the A- weighted loudness to varying flight conditions and some propagation 
parameters. It is to be noted that these plots are dependent on the underlying off-body pressure waveform, 
and hence dependent on the concept under investigation. Some of the trends may be different for other 
near-field waveforms. Several observations may be made from these plots. 

1. The sensitivities with respect to Mach number and step size (A a) are much larger than for cruise 
altitude and number of points used in propagation. Sensitivity with respect to Mach number decreases 
as Mach number increases. This is perhaps because as the Mach number increases, the ground signature 
tends to an N-wave. Put another way, the sensitivity is greater when the ground signature is shaped 
and decreases as the signature tends towards N-waves. 

2. According to the US standard atmosphere, 22 a constant lapse rate exists below 36,089 feet, which 
causes the Blokhintzev invariant to vary as the cruise altitude changes. Above this altitude, in the 
tropopause, cruise altitude variations do not change the ambient density or the Blokhintzev invariant. 
Thus, above 36,089 feet, the sensitivity becomes essentially zero. 

3. As the number of points used in propagation increases, the loudness sensitivity asymptotically goes to 
zero, the loudness values converge while the computational time increases. This is a parameter where 
the sensitivity and uncertainty go down simultaneously. 

4. For small step sizes, sensitivity levels are high. As the step size increases, a point is reached where the 
propagation algorithm determines that the user provided step size is too large to avoid multi-valued 
pressure profiles, and hence picks the largest step size that allows the propagation to run. At this 
step size, the sensitivity value goes to zero. However, the error in loudness might be large because as 
the step size increases, the error induced also increases at each step and the cumulative effect during 
propagation to the ground can be significant. 
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Figure 5: Typical convergence history of the gradient optimizer. 



Figure 6: Ground signatures corresponding to targets. 


Figure 8 shows the sensitivity of the scaling factor k n to Mach number at different conditions. The 
primary take-away from this plot is that the scaling factor sensitivity, and consequently the ray tube area 
sensitivity is small at smaller step sizes compared to larger step sizes. 

This paper considered only some sources of epistemic uncertainty. One possible extension of this work is 
to include sources of aleatory uncertainty such as wind, temperature and relative humidity profiles during 
atmospheric propagation. The foregoing analysis and its associated sensitivities may be used to quantify 
uncertainty during the propagation process. Uncertainty analysis could either be sampling-based, such as 
Monte-Carlo analysis, or non-sampling-based, such as polynomial chaos. Whatever method is chosen, the 
design intent is to reach a robust point where the error as well as sensitivity are simultaneously minimized. 
Even if such a point cannot be reached, the information obtained could be valuable in the understanding 
and exploration of the design space. 
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Figure 7: Sensitivities of flight conditions and selected propagation parameters. 



Figure 8: Sensitivities of scaling factor to Mach number. 
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V. Conclusions 


The adjoint formulation within the sBOOM propagation has been extended to generate sensitivities with 
respect to equivalent areas, flight conditions and selected propagation parameters. Equivalent area targets 
were generated using a gradient-based optimizer at multiple azimuthal angles that minimize ground loudness 
at each off-track location. These may be used during shape optimization exercises for matching at under- and 
off-track locations. The sensitivity of the loudness to flight conditions and propagation parameters has been 
plotted and observations have been made. As the design evolves into the preliminary design stage, analysis 
using sensitivities with respect to flight conditions, propagation parameters and atmospheric properties may 
be used to understand the uncertainty of the underlying propagation analysis as well as provide feedback to 
steer towards a robust design. 


Appendix 

The tridiagonal matrices for the relaxation processes are: 


,A% = 


(l 0 

0 1 0 

0 —OLK\—K2 (1 + 2olKi) K 2 — QLK,\ 


B n ,B% = 


1 0 

0 1 0 

OL , K\—K2 (1 — 2a! K\) + Ob' K\ 


0 

1 / 


0 1 0 

0 1 J 


In the above matrices, , k 2 = Af > and a' = 1 — a. If using the Crank-Nicholson scheme, 

a = 0.5. For thermo-viscous absorption, the matrices are given below with A = 2 r(A r) 2 
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